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Abstract 

Corrections of order a^{Zaf'm to the Lamb shift and corrections of order 
a^[Za)Ep to hyperfine sphtting generated by the insertions of the three- loop 
one-particle reducible diagrams with radiative photons in the electron line are 
calculated. The calculations are performed in the Yennie gauge. 

I. INTRODUCTION 

The theory of Lamb shift in hght hydrogenhke atoms is rapidly developing. Calcula- 
tions of the last corrections which are significantly larger than 1 kHz for the IS* state in 
hydrogen were completed recently. These are corrections of orders a^{Za)'^m, a{Za)"m, 
and a'^{Za)^m. The previously unknown correction to the Lamb shift of order a^{Za)'^m 
is connected with the three-loop contribution to the slope of the Dirac form factor, and 
was obtained in [1]. Corrections of orders a{Za)""m for n = 4,5,6 are already well known 
perturbatively for some time (see, e.g., review [2]), but relatively large magnitude of the 
contributions with n = 6, and high precision of the experimental data [3,4] required calcula- 
tion of the corrections of higher order in Za. All such corrections were obtained numerically 
without expansion in Za [5,6]. Corrections of order a{Za)^m for higher energy levels were 
also calculated recently with high accuracy [7,8]. From the practical point of view the cor- 
rections of order a{Za)^m are no more a significant source of theoretical uncertainty for the 
Lamb shift, for example for the lowest S and P states in hydrogen these contributions are 
now known with uncertainty about 1 Hz. 
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Another recent success is connected with the corrections of order a'^(Za)^m. The leading 
logarithm cubed correction of this order was calculated in the pioneering work [9] . This paper 
was followed by a heated discussion about the magnitude of the nonleading logarithmic 
contributions of order a'^{Za)^m, and even the magnitude of the leading logarithm cubed 
term was put under suspicion (see, e.g., review [2]). The doubts about the magnitude of the 
leading logarithm cubed term were put to rest in [10], where all logarithmically enhanced 
terms of order a'^{Za)^'m were calculated in the Coulomb gauge. Finally, the dominant part 
of the nonlogarithmic contribution of order a'^[Za)^m was obtained in [11]. At the present 
stage remaining uncertainty of the contributions of order Q;^(ZQ;)^m to the Lamb shift is 
about 0.9 kHz and 0.1 kHz for the 15" and 25" states in hydrogen, respectively. 

The largest still unknown contributions to the Lamb shift in hydrogen are corrections of 
order a^{Za)^m. The magnitude of these corrections can be easily estimated multiplying 
the corrections of order a'^{Za)^m [14-16] by an extra factor a/vr. As a result of this simple 
exercise we see that the corrections of order a^{Za)^m should be about 1 kHz for the l-S" 
state in hydrogen. In [17] we calculated radiative corrections to the Lamb shift of order 
a^{Za)^m and radiative corrections to hyperfine splitting of order a^{Za)EF generated by 
the diagrams with insertions of radiative photons and electron polarization loops in the 
graphs with two external photons. 

Below we calculate in the Yennie gauge corrections of order a^{Za)^'m to the Lamb shift 
and corrections of order a^{Za)Ep to hyperfine sphtting generated by the diagrams in Fig. 
1 with insertions of the three-loop one-particle reducible diagrams with radiative photons in 
the electron line. 
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FIG. 1. Reducible three- loop diagrams 



II. FACTORIZED CONTRIBUTIONS TO THE LAMB SHIFT AND HYPERFINE 

SPLITTING 

A. Skeleton Diagram Contributions 

The diagrams for nonrecoil corrections of order a^(Za)^m to the Lamb shift and cor- 
rections of order a^{Za)EF to hyperfine splitting in Fig. 1 can be obtained by three- loop 
radiative insertions in the skeleton diagram in Fig. 2. Respective corrections of lower orders 
in a generated by one- and two-loop radiative insertions are already well known (see, e.g., 
review [2]). All corrections of order a^iZotfm and a^{Za)Ep may be calculated in the 
scattering approximation (see, e.g., [13]). 

Calculation of all these contributions starts with the skeleton diagram in Fig. 2. Contri- 
bution of each of the diagrams in Fig. 1 to the Lamb shift is described by the integral 
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§— — m -—L{k)6io, 1 

Trw^ \m J Jo |kp 

where m is the electron mass, M is the proton mass, = m/(l + m/M) is the reduced 
mass, a is the fine structure constant, Z is the nucleus charge in terms of the electron 
charge {Z = 1 for hydrogen and muonium), and |k| is the magnitude of the dimensionless 
momentum of the external photons measured in the units of the electron mass. The function 
L(k) describes radiative corrections to the skeleton diagram, and should be calculated for 
each particular diagram in Fig. 1. It is normalized to the contribution of the skeleton 
numerator, Lsfee/(k) = 1. The skeleton contribution to the Lamb shift is infrared divergent. 
For some diagrams in Fig. 1 this infrared divergence survives in the Feynman gauge even 
after insertion of the factor L(k) which describes radiative insertions. In order to avoid such 
spurious infrared divergencies which anyway cancel in the gauge-invariant sets of diagrams 
we use infrared safe Yennie gauge for radiative photons in the calculations below. 

Contribution of each of the diagrams in Fig. 1 to hyperfine splitting is described by the 
integral^ 

Trn'^ ^0 |kp 

where -F(k) describes radiative corrections to the skeleton diagram, and, as the function 
L(k) in the case of Lamb shift, should be calculated for each particular diagram in Fig. 1. 
It is normalized to the contribution of the skeleton numerator, Fskei{^) — 1- 



FIG. 2. Skeleton two-photon diagram 



B. Mass Operator and One-Loop Vertex with One On-Mass-Shell Leg in the Yennie 

Gauge 

To calculate contributions of the diagrams in Fig. 1 we use explicit expressions for the 
electron mass operator and the electron-photon vertex in the Yennie gauge. The one-loop 



We define the Fermi energy Ep as 

^F = f^V^(l + a,)(^)%/.i?oo, (2) 

where m is the electron mass, M is the muon mass, nir = m/{l + m/M) is the reduced mass, c is 
the velocity of light, h is the Planck constant, is the Rydberg constant, and is the muon 
anomalous magnetic moment. 
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electron self-energy in the Yennie gauge renormalized on the mass-shell is well known, and 
has the form (see, e.g. [13]) 



(4) 



where 



M(k) 



l-k2 fl-k2)2 



Ink^ 



(5) 



all momenta are dimcnsionless, and are measured in the units of the electron mass m — 1, 
and = (1, 0),pk = 0, k^ = (0, k). 

Renormalized vertex operator with one on-mass-shell leg in the Yennie gauge has the 
form [18] (we omit terms proportional to the momentum k^ because they do not contribute 
to the Lamb shift and hyperfine splitting) 



= ^{A{k)k^-f^ + 5(k)7^(p - A; - 1) + C{k)p^{p - k - 1) + ^(k)a^,A;^}, (6) 



where 



/2 1 \ ^„ ^ 2 ^„ , 3,,,,, ^Ink^ 3 lnk2 



(7) 



C(k) = ^<|.(k) - S{k) - iM(k) + ilnk^ - ^i^. 



(9) 



and 



E{k) = -M$(k) - i5(k) - iM(k) + J + ^Ink^ 



1 d;^ , l+k2^(l-;2) 

m — 



«(k)^ik|/^ 



k^z 



3 lnk2 
81 -k2' 



(10) 



(11) 



Li(l-|k|)-Li(l + |k|) + 2 



Li 1 + 



Vk2 + 1 + Ikl 



-Li 1 



Vk2 + 1 + Ikl 



TT 

T 



^/Wn, vl?+4+|k| 
•J k = — — — — m , = — —. 
^ ' 2|k| Vk2T4-|k| 

Euler dilogarithm Li is defined here as in [13], and the function $(k) usually arises in 
calculations of the diagrams with factorized radiative insertions in the electron line, see, e.g., 
[19]. 
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C. Factorized Corrections of Order o?(ZoLfm, to the Lamb Shift 

Wc use explicit expression for the self-energy operator in eq.(4), and calculating the 
spinor projection on the Lamb shift obtain the contribution of the graph a in Fig. 1 to the 
Lamb shift^ 



1024 647 TT^n^ 
Performing similar calculations for the graph h in Fig. 1 we obtain 



153 , 63\ o?(Zaf /mA^ 



512 327 7r2n3 Vm 
For the diagram c in Fig. 1 we obtain 

a^^Zaf /mr\^ ( 3 



(14) 



+ [B{k) + C{k) - E{k)f\ = -4. 305 82 (1) ( Ih.] m. 



D. Factorized Corrections of Order a^{Za)Ef to Hyperfine Splitting 

Like in the case of the Lamb shift we use explicit expression for the self-energy operator 
in eq.(4), and calculating the spinor projection on hyperfine splitting obtain the contribution 
of the graph a in Fig. 1 to hyperfine splitting 

^^HFS - E^^^ f dk{l - k') (2 - k') M\k) (15) 

/1215 2 81 \a^{ZQ) ^ 
:TT — — I — „ ^ Ep. 



\102A " 64 J n^n^ 



^From now on Ikl = k. 
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Performing similar calculations for the graph b in Fig. 1 we obtain 



a^{Za) 9 



TT^n^ 27r^ Jo 



oc 



2AE''hfs = ' / Ep— / dkM\k) (1 - A;') - A;M(A;) + 2B{k) + C(A;) - E{k) (16) 



, , a^(Za) 
-22. 064 414 1 ; / Ef. 



For the graph c in Fig. 1 we obtain 



^E'hfs = ^ [ ? A dkM(k)lk^A(k) \k^A{k) - 2B{k) - C{k)\ (17) 
+ (2 - A;2) 5(A;) [5(A;) + C{k) - E{k)] + E{k) [B{k) + C(fc) - E{k)] 



, , a^(Za) 
= 11. 723 748 6 ; / 



III. SUMMARY 



In this paper we calculated factorized corrections of order a^{Za)^m to the Lamb shift, 
and factorized corrections of order a^(Za)EF to hyperfine splitting generated by the dia- 
grams in Fig. 1. Collecting contributions to the Lamb shift in eq.(12), eq.(13), and eq.(14), 
we obtain 



AEr = -5. 321 93 (1) 
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m, (18) 



or 

?tot 



AE'f = -0. 535 kHz (19) 



for the IS level in hydrogen. 

Collecting all contributions to hyperfine splitting in eq.(15), eq.(16), and eq.(17) we 
obtain 

~2 



AE'^'fs = 0. 104 23 (1) Ef, (20) 



TT 

or 

?tot 



^Ehfs = 0- 000 13 kHz (21) 

for the ground state in muonium. 

The result in eq.(19) has just the scale we expected on the basis of the general consider- 
ations explained in the Introduction, and corrections of this magnitude are phenomenologi- 
cally relevant at the current level of the experimental and theoretical accuracy (see, e.g. [2]). 
Work on calculation of nonfactorizable contributions is now in progress, and we postpone 
discussion of the phenomenological implications of the results in eq.(19) and eq.(20) until 
its completion. 
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